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One of the most commonly noted 
uncertainties with respect to the modelling of 
cirrus clouds and their effect upon the 
planetary radiation balance is the disputed 
validity of the use of Mie scattering results as 
an approximation to the scattering results of 
the hexagonal plates and columns found in 
cirrus clouds. This approximation has 
historically been a kind of default, a result of 
the lack of an appropriate analytical solution of 
Maxwell’s equations to particles other than 
infinite cylinders and spheroids. Recently, 
however, the use of such approximate 
techniques as the Discrete Dipole 
Approximation has made scattering solutions 
on such particles a computationally intensive 
but feasible possibility. In this study, the 
Discrete Dipole Approximation (DDA) 
developed by Flatau (1992) is used to find such 



Figure 1. Lines of constant a eff . Contour interval is .5 
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solutions for homogeneous, circular cylinders and disks. This 
can serve to not only assess the 
validity of the current radiative 
transfer schemes which are available 
for the study of cirrus but also to 
extend the current approximation of 
equivalent spheres to an 
approximation of second order, 
homogeneous finite circular cylinders 
and disks. The results will be 
presented in the form of a single 
"variable, the extinction efficiency, Q cxt - 
Before proceeding, a few 
definitions and distinctions must be 
given. The first concents the 
normalization of the extinction 
efficiency for nonspherical particles. 
Usually, Q exl is normalized according 
to the following equation: 


Figure Z Ratio of Renormalized Q ext values to Q ext values of corresponding 
infinite cylinders. The contour interval is .1, m * 1.32 + Oi, and 0*90. 

where C ext is the extinction cross section and a eff is the (effective) 
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radius of the particle. This 
method is, of course, 
inappropriate for nonspherical 
particles and in the case of an 
infinite cylinder becomes 
inapplicable. For the infinite 
cylinder case, as defined in 
Bohren and Huffman (1983), 
the normalization becomes: 



where r^ is the cylinder 
radius and L is the cylinder 
length. Of course, this 
represents a normalization 
with respect to actual cross 
sectional area (with 
perpendicular incidence) 
instead of an "effective" cross 
sectional area. 

All of the graphs 
shown in the analysis section 
are contour plots in which 
results are expressed in terms 
of cylinder radius and aspect 
ratio, /3, where: 



It should also be noted that in 
all cases the wavelength of the 
incident radiation is taken to 
be 27 t and that results will be 
given for two common 
refractive indices of ice. 
Since: 
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and all the axes are logarithmic, lines of constant a eff will be linear with a slope of -1/3 (see Figure 1). Also, 
the top and bottom lines in Figure 1 reveal the region of the available data. These appear in all of the plots. 
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2.0 Analysis 

It is most appropriate 
to begin the analysis with a 
look at the Q extl values. With 
the aforementioned 
renormalization we will be 
able to address the degree to 
which finite cylinders behave 
like infinite cylinders. Figures 
2 through 4 display ratios of 
extinction efficiency for finite 
cylinders and infinite cylinders 
of the same radius. Here m, 
equals 1.32 + Oi and the 
angle of incidence (0) is 
expressed with respect to the 
cylinder axis. Note first that a 
"plateau" of infinite cylinder 
behavior is visible in Figure 2 
in a region where r cyi is 
between .4 and 3 and where 
is between about .6 and 10, 
revealing that, in some cases 


t a = 33 . 5573 


at this incidence angle, the behavior of finite cylinders can approach that of their infinite counterparts very 
quickly, sometimes with aspect ratios as low as one or less. The region of this approach shows a strong but 
apparently monotonic dependence on r ,, appearing at lower values of & as r^ increases. There is also a 
dependence upon incidence angle, with the oblique incidence case showing a more complex approach to 
infinitum. Along these lines, resonance peaks are present near ^=2, revealing most likely a Struve-function- 
type approach to infinitum 1 
instead of a purely monotonic 
approach. The peak found in 
the region of ^ = 10 is very 
strong and is most strongly 
observed as 6 approaches 0. 

More on this in a moment. 

Extending this to a case of 
strong absorption (m 3 = 1.32 
+ 0.05 i), the same basic 
structure is observed in Figure 
4, where it becomes apparent 
that the approach to infinitum 
and its dependence upon r^, 
is in fact not monotonic, but 
that in regions below r^, = 

1.0 a strong dependence on 
r^ starts to become visible. A 
quick glance at figure 7 
reveals this to be a region of 
mostly "Mie-type" behavior, 
which in the realm of 
extinction-per-unit-length 
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Figure 6. Same as figure 5 but for ©=33.5573. 
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apparently translates to a high dependence on r , and a relatively low dependence on cylinder length. 

Figures 5 through 7 show ratios of traditionally normalized Q exl values with respect to those of 
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equivalent spheres. Here, in 
the same region around 
^=10, efficiencies which are 
much less than that of spheres 
are observed. A quick glance 
at figure 1 reveals that the 
Tpyj* 10 resonance displayed in 
figures 2 through 4 occurs in 
a region in which a eff is close 
to 4, where the Mie extinction 
profiles reach a maximum. 

Apparently the strong 
resonances in figures 2 
through 4 are due to the disks 
behaving to some extent as 
Mie particles, behavior which 
translates into a large amount 
of extinction per unit length 
in the case of thin disks. In 
contrast, however, the peaks 
near ^,=2 are reiterated 
here, revealing that neither 

spherical nor infinite cylinder 

theory can explain their existence, and that they are most likely resonances associated with characteristic 
distances in the particle. What is most observable is a large degree of remarkable agreement between finite 
cylinders and spheres for a large range of aspect ratios, cylinder radii, incidence angles, and refractive 
indices. As would be expected, the region of agreement is centered around 0=1 and improves slightly with 
the addition of absorption. 

3.0 Conclusion 

The agreement of extinction efficiencies between finite cylinders and spheres is remarkable especially 
at aspect ratios close to 1. This agreement subsides as # increases and the particle begins to behave like 
an infinite cylinder. This agreement also subsides as 0 becomes much less than 1. This agreement region 
shows a high sensitivity to incidence angle and a weak dependence on absorption. Infinite cylinder behavior 
is observed for nonabsorbing cylinders at normal incidence in some cases even when /? is less than 1, but 
this agreement becomes more complex as incidence angles become oblique. 
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Figure 7. Same as figure 6, but for m=“1.32 + .05i. 
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